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Recent observation of high-pr hadron spectra suppression and mono-jet production in central 
Au + Au collisions and their absence in d + Au collisions at RHIC have confirmed the long pre- 
dicted phenomenon of jet quenching in high-energy heavy-ion collisions. Detailed analyses of the 
experimental data also show parton energy loss as the mechanism for the discovered jet quenching. 
Using a pQCD parton model that incorporates medium modified parton fragmentation functions 
and comparing to experimental data from deeply inelastic scattering off nuclei, one can conclude 
that the initial gluon (energy) density of the hot matter produced in central Au + Au collisions that 
causes jet quenching at RHIC is about 30 (100) times higher than in a cold Au nucleus. Combined 
with data on bulk and collective properties of the hot matter, the observed jet quenching provides 
strong evidence for the formation of a strongly interacting quark gluon plasma in central Au + Au 
collisions at RHIC. 



I. INTRODUCTION 

In the last three years, since the physics opera- 
tion of the Relativistic Heavy-ion Collider (RHIC) at 
Brookhaven National Laboratory (BNL), one has wit- 
nessed tremendous progress in heavy-ion experimental 
physics and many milestones in the search for the elusive 
quark gluon plasma (QGP) that is expected to be formed 
in high-energy heavy-ion collisions. One of the most ex- 
citing phenomena observed at RHIC is jet quenching, 
or suppression of leading hadrons from fragmentation of 
hard partons due to their strong interaction with the 
dense medium. Such a phenomenon was long predicted 
by a pQCD-based model calculation [1] , but was not seen 
in high-energy heavy-ion collisions until the first mea- 
surements in central Au + Au collisions at yfs = 130 GeV 
at RHIC [2,3]. The azimuthal distribution of high pt 
hadrons was also found to display large anisotropy with 
respect to the reaction planes [4] of non-central Au + Au 
collisions which was expected to be caused by the path- 
length dependence of jet quenching [5,6]. 

The measurements were confirmed and further ex- 
tended to a larger px range in Au + Au collisions at 
y/s = 200 GeV, and the suppression was found to have 
a weak p-r-dependence at pt > 6 GeV/c [7,8], indepen- 
dent of hadron species. At the same time, the back- 
side high-j»T two-hadron correlation in azimuthal angle, 
characteristic of high-py back-to-back jets in p + p colli- 
sions, was found to disappear in central Au+Au collisions 
[9], confirming the predicted mono-jet phenomenon of jet 
quenching [10]. 

The final milestone in the experimental discovery of 
jet quenching was achieved during the third year of the 
RHIC physics program, when both the single-hadron 
spectra [11-13] and the disappearance of away-side two- 
hadron correlation [12] were found to be absent in the 
same central rapidity region of d + Au collisions at 
y/s = 200 GeV. These d + Au results prove that the ob- 
served high-pT suppression patterns in Au+Au collisions 



are not initial state effects encoded in the wavefunction 
of a beam nucleus, but are jet quenching caused by fi- 
nal state interaction of hard partons with the produced 
dense medium. 

Following these three major experimental results, ad- 
ditional data — such as the dependence of away-side sup- 
pression on the azimuthal angle relative to the reaction 
plane [14], modification of the hadron distributions (frag- 
mentation functions) both along and opposite the direc- 
tion of the triggered high-p^ hadron [15], and absence of 
suppression in the direct photon spectra [16] in Au + Au 
collisions — have now solidified the conclusion that the 
observed jet quenching is caused by parton multiple scat- 
tering and energy loss in the hot and dense medium. 
Furthermore, the splitting of baryon and meson spec- 
tra suppression and azimuthal anisotropy in the interme- 
diate pt < 6 GeV/c region [17-19] also point to non- 
trivial medium modification of hadronization that is in- 
dicative of dense partonic matter. Such a wealth of data 
affords one a quantitative phenomenological analysis of 
jet quenching and a tomographical picture of the hot and 
dense matter formed in heavy-ion collisions at RHIC. 

In this report, we will first briefly review the recent 
progress in pQCD study of parton multiple scattering 
and induced radiative energy loss in a dense medium. 
Combining this with the pQCD parton model of high- 
Pt jet and hadron production, one can analyze the ob- 
served jet quenching phenomena to extract properties of 
the dense matter produced. One can combine the de- 
duced properties with other information from analyses 
of bulk particle spectra such as collective flow and total 
multiplicity and energy production to present a collection 
of compelling evidence for the existence of a strongly in- 
teracting quark gluon plasma in Au + Au collisions at 
RHIC. However, we will focus in this report only on the 
analysis of jet quenching data and the properties of the 
dense matter one can extract from these data. 
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II. PARTON ENERGY LOSS 

One important step in establishing evidence of QGP 
formation is to characterize the properties of the pro- 
duced dense medium, for example, the parton and en- 
ergy density and current-current correlations, among 
many other characteristics. Traditionally, one can study 
the properties of a medium via scattering experiments 
with particles. In deeply inelastic scattering (DIS) ex- 
periments, for example, leptons scatter off the nucleon 
medium via photon exchange. The response function or 
the correlation function of the electromagnetic currents 

W^(q) = ^J d 4 xe^ x (A | j; m (0)jt m (x) I A), (1) 

is a direct measurement of the quark distributions in a 
nucleon or nucleus. For a dynamic system in heavy-ion 
collisions, one can no longer use the technique of scatter- 
ing with an external beam of particles in the conventional 
sense because of the transient nature of the matter. The 
lifetime of the system is very short, on the order of a 
few fm/c. The initial spatial size is only the size of the 
colliding nuclei, about 6 fm in diameter in the transverse 
dimension for the heaviest nuclei. The system expands 
very rapidly both in the longitudinal and the transverse 
direction. These characteristics make it impossible to use 
external probes to study the properties of the produced 
dense matter in high-energy heavy- ion collisions. One 
has to resort to internal probes such as the electromag- 
netic emission from the dynamic medium, whose rate is 
given by the thermal average of the above correlation 
function. The emission rate depends mainly on the local 
temperature or the parton density while the total yield 
also depends on the whole evolution history of the sys- 
tem. Therefore, a strongly interacting system can reveal 
its properties and dynamics through photon and dilep- 
ton emission. The same medium interaction should also 
cause attenuation of fast and energetic partons propagat- 
ing through the medium. Such an effect is the underlying 
physics of the jet quenching [20,21] phenomenon and jet 
tomography technique for studying properties of dense 
matter in high-energy heavy-ion collisions. 

Jet quenching as a probe of the dense matter in heavy- 
ion collisions takes advantage of the hard processes of jet 
production in high-energy heavy-ion collisions. Because 
large-j»T partons are produced very early in heavy-ion 
collisions, they can probe the early stage of the formed 
dense medium. Similar to the technology of computed 
tomography (CT) , the study of these energetic particles, 
their initial production and interaction with the dense 
medium can yield critical information about the proper- 
ties of the medium that is otherwise difficult to access 
through soft hadrons from the hadronization of the bulk 
matter. Though relatively rare with small cross sections, 
the jet production rate can be calculated perturbatively 
in QCD and agrees well with experimental measurements 



in high-energy p + p(p) collisions. A critical component 
of the jet tomography is then to understand the jet's at- 
tenuation through dense matter as it propagates through 
the medium. 

There have been many theoretical studies [22-24] of jet 
attenuation in a hot medium in recent years. The first 
attempt was by Bjorken [20] to calculate elastic energy 
loss of a parton via elastic scattering in the hot medium. 
A simple estimate can be given by the thermal averaged 
energy transfer v c \ « q^/2uj of the jet parton to a ther- 
mal parton, with energy u, and q± the transverse mo- 
mentum transfer of the elastic scattering. The resultant 
clastic energy loss [25] 

is sensitive to the temperature of the thermal medium but 
is small compared to radiative energy loss. Here, \i is the 
Debye screening mass and C2 is the Casimir coefficient of 
the propagating parton in its fundamental presentation. 
The elastic energy loss can also be calculated within fi- 
nite temperature field theory of QCD [26] with a more 
careful and consistent treatment of the Debye screening 
effect. 

Though there had been estimates of the radiative par- 
ton energy loss using the uncertainty principle [27], the 
first theoretical study of QCD radiative parton energy 
loss incorporating Landau- Pomeranchuk-Migdal interfer- 
ence effect [28] is by Gyulassy and Wang [29] where mul- 
tiple parton scattering is modeled by a screened Coulomb 
potential model. Baier et al. (BDMPS) [30] later consid- 
ered the effect of gluon rescattering which turned out to 
be very important for gluon radiation induced by multi- 
ple scattering in a dense medium. These two studies have 
ushered in many recent works on the subject, including a 
path integral approach to the problem [31], twist [32] or 
opacity expansion framework [33,34] which is more suit- 
able for multiple parton scattering in a finite system. The 
radiative parton energy loss in the leading order of the 
twist or opacity expansion was found to have a simple 
form [33] 

_« Ca — -]!!(-_), (3) 

for an energetic parton (E > (i) in a static medium, 
where A is the gluon's mean free path in the medium. The 
unique L-dcpendcnce of the parton energy loss, which 
was first discovered by BDMPS [30] is a consequence 
of the non-Abclian LMP interference effect in a QCD 
medium. In a dynamic system the total energy loss was 
found [6,35,36] to be proportional to a line integral of the 
gluon density along the parton propagation path. 

Since gluons are bosons, there should also be stim- 
ulated gluon emission and absorption by the propagat- 
ing parton in the presence of thermal gluons in the hot 
medium. Such detailed balance is crucial for parton ther- 
malization and should also be important for calculating 
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the energy loss of an energetic parton in a hot medium. 
Taking into account such detailed balance in gluon emis- 
sion, Wang and Wang [37] obtained an asymptotic behav- 
ior of the effective energy loss in the opacity expansion 
framework , 



AE 
~E 



AXE 
t s C F LT 2 
~3 XE 2 



In 



In 



2E 



0.048 



f 



- 7E 



6C'(2) 



(4) 



where the first term is from the induced bremsstrahlung 
and the second is due to gluon absorption in detailed 
balance, which effectively reduces the total parton en- 
ergy loss in the medium. Though the effect of detailed 
balance is small at large energy (E ^> fi), it changes the 
effective energy dependence of the energy loss in the in- 
termediate energy region. 



III. MODIFIED FRAGMENTATION FUNCTION 

Because of color confinement in the vacuum, one can 
never separate hadrons fragmenting from the leading par- 
ton and particles materializing from the radiated gluons. 
The total energy in the conventionally defined jet cone 
in principle should not change due to induced radiation, 
assuming that most of the energy carried by radiative glu- 
ons remains inside the jet cone [38]. Additional rescatter- 
ing of the emitted gluon with the medium could broaden 
the jet cone significantly, thus reducing the energy in a 
fixed cone. However, fluctuation of the underlying back- 
ground in high-energy heavy-ion collisions makes it very 
difficult, if not impossible, to determine the energy of a 
jet on an event-by-event base with sufficient precision to 
discern a finite energy loss of the order of 10 GeV. Since 
high-pT hadrons in hadronic and nuclear collisions come 
from fragmentation of high-pT jets, energy loss naturally 
leads to suppression of high-py hadron spectra. This was 
why Gyulassy and Wang proposed [1] to measure the sup- 
pression of high-pT hadrons to study parton energy loss 
in heavy-ion collisions. 

Since parton energy loss effectively slows down the 
leading parton in a jet, a direct manifestation of jet 
quenching is the modification of the jet fragmentation 
function, D a —yh(z, /jt), which can be measured directly 
in events in which one can identify the jet via a com- 
panion particle like a direct photon [39] or a triggered 
high pt hadron. This modification can be directly trans- 
lated into the energy loss of the leading parton. Since 
inclusive hadron spectra are a convolution of the jet pro- 
duction cross section and the jet fragmentation function 
in pQCD, the suppression of inclusive high-px hadron 
spectra is a direct consequence of the medium modifica- 
tion of the jet fragmentation function caused by parton 
energy loss. 

Since a jet parton is always produced via a hard pro- 
cess involving a large momentum scale, it should also 



have final state radiation with and without rescatter- 
ing, leading to the DGLAP evolution of fragmentation 
functions. Such final state radiation effectively acts as a 
self-quenching mechanism, softening the leading hadron 
distribution. This process is quite similar to the induced 
gluon radiation and the two should have a strong interfer- 
ence effect. It is therefore natural to study jet quenching 
and modified fragmentation functions in the framework 
of modified DGLAP evolution equations in a medium 
[32]. 

To demonstrate medium modified fragmentation func- 
tion and parton energy loss, we review here the twist 
expansion approach in the study of deeply inelastic scat- 
tering (DIS) eA [32,40]. While the results can be readily 
extended to the case of a hot medium, they also provide 
a baseline for comparison to parton energy loss in cold 
nuclei. 

In the parton model with the collincar factorization ap- 
proximation, the leading-twist contribution to the semi- 
inclusive cross section can be factorized into a product of 
parton distributions, parton fragmentation functions and 
the hard partonic cross section. Including all leading log 
radiative corrections, the lowest order contribution from 
a single hard 7* + q scattering can be written as 



dz h 



^2 e l J dxfq{x,fj, 2 I )H^J(x,p,q)D q ^ h (z h ,fM 2 ), 



where HjxJ (x,p,q) is the hard part of the process in 
leading order, £h the observed hadron momentum, p = 
[j> + ,0,Oj_] the momentum per nucleon in the nucleus, 
q = [— Q 2 /2q~, q~ , 0±] the momentum transfer carried 
by the virtual photon. In the chosen frame, q~ is the 
quark momentum transferred from the virtual photon. 
The momentum fraction carried by the hadron is defined 
as Zh = ^h/q~ and x = xb = Q 2 /2p + q~ is the Bjorken 
variable, /i 2 and /j, 2 are the factorization scales for the 
initial quark distributions f£(x, /J, 2 ) in a nucleus and the 
fragmentation functions D q ^} l {zh 1 yu 2 ), respectively. 

In a nuclear medium, the propagating quark in DIS 
will experience additional scattering with other partons 
from the nucleus. The rescatterings may induce addi- 
tional gluon radiation and cause the leading quark to 
lose energy. Such induced gluon radiations will effectively 
give rise to additional terms in a DGLAP-likc evolution 
equation leading to the modification of the fragmentation 
functions in a medium. These are the so-called higher- 
twist corrections since they involve higher-twist parton 
matrix elements and are power-suppressed. The leading 
contributions involve two-parton correlations from two 
different nucleons inside the nucleus. 

One can apply the generalized factorization [41] to 
these multiple scattering processes. In this approxima- 
tion, the double scattering contribution to radiative cor- 
rection can be calculated and the effective modified frag- 
mentation function is [32] 
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[Aj q ^ qg (z, x, xl, £t) D q ^h(z h /z) 
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in much the same form as the DGLAP correction in 
vacuum, where D q ^h{zh, y?) and D g ^h(zh, ^ 2 ) are the 
leading-twist quark and gluon fragmentation functions. 
The modified splitting functions are given as 
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(6) 
(7) 



Here, the fractional momentum x — xb is carried by the 
initial quark and xl — £\ /2p + q~ z{l — z) is the additional 
momentum fraction carried cither by the quark or gluon 
in the secondary scattering that is required for induced 
gluon radiation. The twist-four parton matrix element of 
the nucleus, 



^— dy-dy-e l{x+XL)p+v (1 
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x(l - e ~ iXLP+ ( y -»i ] )e{-y^ )6{y- - t/f) 



x-(A|^(0) 7 + F (T + (j/ 2 -)F +ff (j / r)V,(j/-)|A) , (8) 

has a dipole-like structure which is a result of LPM in- 
terference in gluon bremsstrahlung. 
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FIG. 1. Predicted nuclear modification of the jet fragmen- 
tation functions compared to the HERMES data [43] on ratios 
of hadron distributions between A and d targets in DIS. 

Averaged over a Gaussian nuclear distribution, the in- 
terference will produce a factor 1 — e~ x ^l x A with xa — 
1/MRa- Here Ra is the nuclear size and M is the 
nucleon mass. Using the factorization approximation 
[32,41,42], one can relate the twist-four parton matrix 



elements of the nucleus to the twist-two parton distribu- 
tions of nucleons and the nucleus, 



T q A g (x,x L )*—(l-e- x l/ x -)f q A (x) 

XA 



(9) 



where C = 2CxTf g N (xT) is considered a constant. One 
can identify 1/xlp + = 2q~ z(\ — z)j£\ as the formation 
time of the emitted gluons. In the limit of collinear radi- 
ation (xl — > 0) or when the formation time of the gluon 
radiation is much larger than the nuclear size, the above 
matrix element vanishes, demonstrating a typical LPM 
interference effect. 

Since the LPM interference suppresses gluon radia- 
tion whose formation time (t/ <~ Q 2 /£ 2 r ,p + ) is larger 
than the nuclear size MRa/p + (here in the infinite mo- 
mentum frame), f, 2 ^ should then have a minimum value 



of 4 



Q 2 /MR A ~ Q 2 /A 1 / 3 . Therefore, the leading 
higher- twist contribution is proportional to (x s RaI £\ ~ 
a s R 2 A /Q 2 due to double scattering and depends quadrat- 
ically on the nuclear size Ra- 

With the assumption of the factorized form of the 
twist-4 nuclear parton matrices, there is only one free 
parameter C(Q 2 ) which represents quark-gluon correla- 
tion strength inside nuclei. Once it is fixed, one can pre- 
dict the z, energy and nuclear dependence of the medium 
modification of the fragmentation function. Shown in 
Fig. 1 are the calculated [35] nuclear modification factor 
of the fragmentation functions for 1A N and s4 Kr targets 
as compared to the recent HERMES data [43]. The pre- 
dicted shape of the z-dependence agrees well with the ex- 
perimental data. A remarkable feature of the prediction 
is the quadratic A 2 / 3 nuclear size dependence, which is 
verified for the first time by an experiment. By fitting the 
overall suppression for one nuclear target, one obtains the 
only parameter in the calculation, C(Q 2 ) = 0.0060 GeV 2 
with a s (Q 2 ) = 0.33 at Q 2 w 3 GeV 2 . The predicted 
energy {v) dependence also agrees with the experimen- 
tal data as shown in Fig. 2. Note that the suppression 
goes away as the quark energy increases. This is in sharp 
contrast to the jet quenching in heavy-ion collisions. 
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FIG. 2. Energy dependence of the nuclear modification 
compared with the HERMES data [43]. 
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One can quantify the modification of the fragmenta- 
tion by the quark energy loss which is defined as the 
momentum fraction carried by the radiated gluon, 
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(10) 



where x^ = p 2 /2p + q~ z(l — z) = xb/z(1 — z) if the fac- 
torization scale is set as p 2 — Q 2 . When xa C € 1 
the leading quark energy loss is approximately 



/a \r 2n p^ c Aa 2 s x b r 1 



^V c Q 2 x 2 A 



2x B 



(11) 



Since xa = 1/MRa, the energy loss (Az g ) thus depends 
quadratically on the nuclear size. 

In the rest frame of the nucleus, p + = M, q~ = 
v, and x B = Q 2 /2p + q~ = Q 2 /2Mv. One can get 
the averaged total energy loss as AE = u(Az g ) rts 
C(Q 2 )a 2 (Q 2 )MR 2 A (C A /N c )3\n(l/2xB). With the deter- 
mined value of C, (x B ) w 0.124 in the HERMES exper- 
iment [43] and the average distance (La) = Ra\/2/tt 
for the assumed Gaussian nuclear distribution, one gets 
the quark energy loss dE/dL rts 0.5 GeV/fm inside a Au 
nucleus. 



IV. JET QUENCHING IN HOT MEDIUM 

To extend the study of modified fragmentation func- 
tions to jets in heavy-ion collisions, one can assume 
(k 2 -) w /j 2 (the Debye screening mass) and a gluon den- 
sity profile p(y) = (to/t)8(Ra — y)po for a 1-dimensional 
expanding system. Since the initial jet production rate 
is independent of the final gluon density, which can be 
related to the parton-gluon scattering cross section [44] 
[a s XTG(xT) ~ M 2(J srL one nas then 



OisTqgiXB^XL) 



p 2 f dya g p(y)[l- cos(y/T f )}, (12) 



where r/ = 2Ez(\ — z)/£^ is the gluon formation time. 
The averaged fractional energy loss is then, 
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Keeping only the dominant contribution and assuming 
o g rts C a 2-Ka 2 / p 2 (C a =l for and 9/4 for 53 scatter- 
ing), one obtains the total energy loss, 



(AE) rts ttC, 
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(14) 



which has also be obtained in the opacity expansion ap- 
proach [6,36] in a thin plasma. 

In a static medium, the gluon density is independent 
of time and one can recover the quadratic lenghth depen- 
dence in Eq. (3). Neglecting the logarithmic dependence 
on r, the averaged energy loss can be expressed as 



dE 

{ dL )ld M ( dE o/d L )(^ /R A ) 



(15) 



in a 1-dimensional expanding system, p(r) = p t /t. 
Here dEo/dL cx pqRa is the energy loss in a static 
medium with the same gluon density p as in the 1-d 
expanding system at time r . Because of the expansion, 
the averaged energy loss (dE/dL) id is suppressed as com- 
pared to the static case and does not depend linearly on 
the system size at a fixed value of initial gluon density. 

One should emphasize now that the above simple log- 
arithmic energy dependence of the parton energy loss in 
Eqs. (11) and (14) is only an asymptotic behavior. How- 
ever, kinematic limits on induced gluon radiation from a 
parton with finite energy can result in a stronger energy 
dependence [33] . The effect of detailed balance with ther- 
mal absorption further increases the energy dependence 
[37], while reducing the effective parton energy loss for 
intermediate values of parton energy. Such a detailed 
balance effect sets the energy dependence of the energy 
loss in cold nuclei in DIS apart from the hot medium in 
heavy-ion collisions. If one parameterizes the energy de- 
pendence of the energy loss including the full kinematic 
limits and thermal absorption, one would get 
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/(7.5 + E/p). 



(16) 



The threshold is the consequence of gluon absorption that 
competes with radiation and effectively shuts off the en- 
ergy loss. The parameter p is set to be 1.5 GeV in the 
following discussions. Such a detailed balance effect can 
also explain why the total hadron multiplicity that is 
dominated by soft hadrons does have significant enhance- 
ment over the p + p collisions as result of induced gluon 
emission. 



A. Single Spectra 

To calculate the modified high-p^ spectra in A + A 
collisions, we use a LO pQCD model [45], 



dyd 2 p T 



= K j d 2 bd 2 rdx a dx b d 2 k a Td 2 k T 

abed 

x t A {r)t A {\h - r\)g A (kaT,r)g A (k b T, |b - r|) 
x fa/A{xa,Q 2 ,r)f b/A (x bl Q 2 ,\b-r\) 
D' h/c (z c ,Q 2 ,AE c ) da 



irz c 



dt 



(ab — > cd), 



(17) 



5 



with medium modified fragmentation functions D' h ,, 
Here, z c = pt/ptc, V = Uc cr(ab — > cd) are elementary 
parton scattering cross sections and tA(b) is the nuclear 
thickness function nor malized to J d 2 bt A (b) = A. The 
K 1.5-2 factor is used to account for higher order 
pQCD corrections. 
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FIG. 3. Single-inclusive spectra of charged hadrons in p + p 
collisions at y/Z = 200, 900, 1800 GeV. The solid (dot-dashed) 
lines are parton model calculations with (without) intrinsic 
kr- Experimental data are from Refs. [50,51]. 

To simplify the incorporation of the medium modifica- 
tion of the fragmentation function in the parton model, 
an effective model [39] can be used: 



D' h/C (z c , Q\ AE C ) = e-^D° h/c (z c , Q 2 ) 



Zc 
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(—)^D° (z' 
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where z' c , z' g are the rescaled momentum fractions after 
parton energy loss. This effective model is found [46] to 
reproduce the pQCD result from Eq.(5) very well, but 
only when Az — AE C /E is set to be Az w 0.6(z g ). 
Therefore the actual averaged parton energy loss should 
be AE/E = 1.6Az, with Az extracted from the effective 
model. The factor 1.6 is mainly caused by the unitarity 
correction effect in the pQCD calculation. The fragmen- 
tation functions in free space D^ c (z c ,Q 2 ) will be given 
by the BBK parameterization [47] . 

The parton distributions per nucleon fa/A( x a,Q 2 ,T) 
inside the nucleus are assumed to be factorizable into 
the parton distributions in a free nucleon given by 



the MRSD— ' parameterization [48] and the impact- 
parameter dependent nuclear modification factor given 
by the new HIJING parameterization [49]. The initial 
transverse momentum distribution gA(kr,Q 2 ,b) is as- 
sumed to have a Gaussian form with a width that in- 
cludes both an intrinsic part in a nucleon and nuclear 
broadening. This parton model can describe high-pT 
hadron spectra in p + p{p) well, as shown in Fig. 3 that 
compares the parton model calculation [45] with the ex- 
perimental data on the inclusive charged hadron spectra 
at different collider energies. 
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FIG. 4. The first prediction [45] of the Cronin effect in 
p + Au collisions at y/s = 200 GeV compared to the recent 
RHIC data from PHENIX [11] and STAR [12]. The dashed 
and dot-dashed lines correspond to HIJING [49] and EKS [52] 
parameterizations of nuclear parton distributions. 

In p + A collisions, this parton model incorporates 
both nuclear modification of the parton distributions and 
broadening of the initial transverse momentum. The 
initial momentum broadening leads to an enhancement, 
known as the Cronin effect, of the hadron spectra in the 
intermediate pt region. The enhancement disappears at 
large px- The parameters have been fitted to the nuclear 
modification of the px spectra in p + A collisions at up 
to the Fermilab energy y/s = 40 GeV. Shown in Fig. 4 
is the first prediction made in 1998 [45] of the Cronin 
effect at RHIC for p + Au collisions at y/s = 200 GeV as 
compared to the RHIC data on d + Au collisions. As one 
can see, the initial multiple scattering in nuclei can give 
some moderate Cronin enhancement of the high-p^ spec- 
tra, though the details depend on the parameterization of 
the nuclear modification of parton distributions. There- 
fore, any observed suppression of the high-pT spectra in 
Au + Au collisions has to be caused by jet quenching. 

In A + A collisions, one has to consider medium 
modification of the fragmentation functions due to par- 
ton energy loss according to Eq. (18). Assuming a 1- 
dimensional expanding medium with a gluon density 
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Pg(r,r) that is proportional to the transverse profile 
of participant nucleons, one can calculate the impact- 
parameter dependence of the energy loss, 



AE(b, 



x dL' 



id 



AL 



dT- 



T 



tqPo 



Pg(r, b,r + nT), 



(19) 



according to Eq. (14), where AL(b,r,(f>) is the distance 
a jet, produced at r, has to travel along n at an az- 
imuthal angle </> relative to the reaction plane in a colli- 
sion with impact-parameter b. Here, po is the averaged 
initial gluon density at initial time tq in a central collision 
and (dE/dL)id is the average parton energy loss over a 
distance R A in a 1-d expanding medium with an initial 
uniform gluon density po . The corresponding energy loss 
in a static medium with a uniform gluon density po over 
a distance R A is [35] dE /dL = (R A /2T Q )(dE/dL) ld . We 
will use the parameterization in Eq. (16) for the effective 
energy dependence of the quark energy loss. 
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FIG. 5. Hadron suppression factors in Au + Au collisions 
as compared to data from STAR [8] and PHENIX [7]. See 
text for a detailed explanation. 

Shown in Fig. 5 are the calculated nuclear modifica- 
tion factors Rab(pt) = da\ B / (N Vl na,vy)d(Jp p for hadron 
spectra (\y\ < 0.5) in Au + Au collisions at ^/s = 200 
GeV, as compared to experimental data [8,7]. Here, 
(^binary) = J d 2 bd 2 rt A (r)t A (\b - r\). To fit the observed 
7T° suppression (solid lines) in the most central collisions, 
we have used p — 1.5 GeV, eo = 1.07 GeV/fm and 
A = l/(crp ) = 0.3 fm in Eqs. (18) and (16). The 
hatched area (also in other figures in this report) indi- 
cates a variation of eq = ±0.3 GeV/fm. The hatched 



boxes around R A b = 1 represent experimental errors of 
STAR data in overall normalization, including a constant 
factor of about 18% from the determination of total p+p 
inelastic cross section. Nuclear kx broadening and parton 
shadowing together give a slight enhancement of hadron 
spectra at intermediate pt = 2 — 4 GeV/c without parton 
energy loss. 




O PHENIX n° p T > 4 GeV 
Au+Au -> tt°, Vs=200 GeV 




50 100 150 200 250 300 350 



<N port > 



FIG. 6. The centrality dependence of the measured single 
inclusive hadron suppression [7,8] at high-pr as compared to 
theoretical calculation with parton energy loss 

The flat px dependence of the 7r° suppression is a con- 
sequence of the strong energy dependence of the parton 
energy loss as also noted in Refs. [53,54]. This is in sharp 
contrast to the energy dependence of the suppression in 
DIS as shown in Fig. 2 and should be a strong indication 
of the detailed balance effect in a thermal medium. Such 
an effect also causes the slight rise of R A b at pt < 4 
GeV/c in the calculation. In this region, one expects the 
fragmentation picture to gradually lose its validity and 
to be taken over by other non-perturbative effects, es- 
pecially for kaons and baryons. As a consequence, the 
(K + p)/w ratio in central Au + Au collisions is signifi- 
cantly larger than in peripheral Au+Au or p+p collisions 
[17]. This has been attributed to parton coalescence in 
the fragmentation of the leading parton [55-57]. The 
slight flavor dependence of the Cronin effect in d + Au 
collisions as shown in Fig. 4 could also be attributed to 
the coalescence effect [59,58]. Such an effect, while pro- 
viding interesting insight into the hadronization mecha- 
nism in the quark gluon plasma, will complicate the pic- 
ture of jet quenching and introduce large uncertainties 
in the physics extracted from jet tomography analysis of 
the experimental data. Fortunately, the effect has been 
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shown to disappear at large px > 5 GeV/c both in the 
coalescence models and the experimental data. The sup- 
pression ratio of charged hadrons and 7r° and of A and K 
all converge [19] at large pr- In the calculation shown in 
Fig. 5, a nuclear-dependent (proportional to (Aginary)) 
soft component was added to kaon and baryon fragmen- 
tation functions to take into account the coalescence ef- 
fect, so that (K + p)/n r* 2 at pr ~ 3 GeV/c in the most 
central Au+Au collisions and it approaches its p+p value 
at pt > 5 GeV/c. To demonstrate the sensitivity to the 
parameterized parton energy loss in the intermediate pt 
region, we also show R AA in 0-5% centrality (dashed line) 
for n = 2.0 GeV and e = 2.04 GcV/fm without the soft 
component. 

The measured centrality dependence of the single 
hadron suppression in Au+Au collisions, shown in Fig. 6, 
agrees very well with parton model with energy loss. This 
is the consequence of the centrality dependence of the 
averaged total energy loss in Eq. (14), which leads to 
an effective surface emission of the surviving jets. Jets 
produced around the core of the overlapped region are 
strongly suppressed, since they lose the largest amount 
of energy. The centrality dependence of the suppression 
is found to be more dominated by the geometry of the 
produced dense matter than the length dependence of 
the parton energy loss. 
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FIG. 7. (a) Two-particle azimuthal distributions for mini- 
mum bias and central d + Au collisions, and for p+p collisions 
[12]. (b) Comparison of two-particle azimuthal distributions 
for central d+Au collisions to those seen in p+p and central 
Au+Au collisions. 



B. Di- hadron Spectra 

Since jets are always produced in pairs in LO pQCD, 
two-hadron correlations in p + p collisions should have 
unique characteristics pertaining to the back-to-back jet 
structure of the initial hard parton-parton scattering. Jet 
quenching should also modify this di-hadron correlation 
of the back-to-back jets. Shown in Fig. 7 are the two- 
hadron correlations in azimuthal angle as measured by 
STAR experiment in p+p, d + Au and Au + Au collisions 



at RHIC [12]. While the same-side correlation remains 
the same, the away-side correlation due to the back-side 
jet is strongly suppressed in central Au + Au collisions. 
The energy loss that causes the suppression of the single 
inclusive hadron spectra should also be able to explain 
the disappearance of away-side correlation. 
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FIG. 8. Back-to-back correlations for charged hadrons with 
p^ iB > p T > 2 GeV/c, p^ ig = 4-6 GeV/c and \y\ < 0.7 in 
Au + Au (lower curves) and p + p (upper curves) collisions as 
compared to the STAR [9] data. 

In the same LO pQCD parton model, one can calculate 
the spectra [60], 



zp jp aa AA 



-E 

2 ^ 

abed 



d bd rdx a dx d k a xd k^T 



d 3 pid 3 p 2 

x t A (r)t A (\b - r\)g A (k aT , r)g A (k bT , |b - r| 
x fa/A(x a ,Q 2 ,r)f b/A {x b ,Q 2 ,\b -r|) 

x D h/c {z c , Q\ AE c )D h/d (z d , Q 2 , AE d ) 

x ^(a6-+ cd)5 4 (p a +p b ~p c -p d ), 
dt 



2nz 2 zj 



(20) 



of two back-to-back hadrons from independent fragmen- 
tation of the back-to-back jets. Let us assume hadron 
h\ is a triggered hadron with pti = p^ Ig . One can de- 
fine a hadron-triggered fragmentation function (FF) as 
the back-to-back correlation with respect to the triggered 
hadron: 

72 tri 



D h ^{z T ^^)=p^ 



d^t/dY^dpTdcji 

similarly to the direct-photon triggered FF [39] in 7- 
jet events. Here, zt — Pt/Pt 1S an d integration over 
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I2/1.2 1 < At/ is implied. In a simple parton model, two 
jets should be exactly back-to-back. The initial parton 
transverse momentum distribution in our model will give 
rise to a Gaussian-like angular distribution. In addition, 
we also take into account transverse momentum smearing 
within a jet using a Gaussian distribution with a width 
of (k±) = 0.6 GeV/c. 

Shown in Fig. 8 are the calculated back-to-back cor- 
relations for charged hadrons in Au + Au collisions as 
compared to STAR data [9]. The same energy loss that 
is used to calculate single hadron suppression can also 
describe well the observed away-side hadron suppression 
and its centrality dependence. 

With cross sections integrated over 0, one obtains 
a hadron-triggered FF D hlfl2 (z T , p^ ie ) . The suppres- 
sion factor I AA {z Tl p t ^) = D^/Dfy 112 defined by the 
STAR experiment [9] is just the medium modification 
factor of the hadron-triggered FF. The shape of the mod- 
ification is predicted [60] to be very similar to that of the 
7-triggered fragmentation functions [39]. 



tra, one will obtain azimuthal anisotropic spectra. Shown 
in Fig. 9 is V2{pt) (second Fourier coefficient of the az- 
imuthal angle distribution) of charged hadrons generated 
from parton energy loss (dot-dashed) as compared to pre- 
liminary STAR data [61], using the 4-particle cumulant 
moments method [62] which is supposed to reduce non- 
geometrical effects such as inherent two-particle corre- 
lations from di-jet production [63]. The energy loss ex- 
tracted from high-pT hadron spectra suppression can also 
account for the observed azimuthal anisotropy at large 
Pt- At intermediate pr, the observed is larger than 
the simple parton model calculation. Such discrepancy 
can also be attributed to effects of parton coalescence. If 
the difference is made up by kaons and baryons from the 
coalescence contribution, one finds that they must have 
v 2 « 0.23 (0.11) for 20-50% (0-10%) collisions. The total 
V2(pt) is shown by the solid lines. 



V. JET QUENCHING AND QGP 



C. High pt Azimuthal Anisotropy 

Another aspect of jet quenching is the azimuthal 
anisotropy of the spectra caused by parton energy loss 
which depends on the azimuthal angle according to 
Eq. (19). In non-central collisions, the average path 
length of parton propagation will vary with the azimuthal 
angle relative to the reaction plane. This leads to an az- 
imuthal dependence of the total parton energy loss and 
therefore azimuthal asymmetry of high-p^ hadron spec- 
tra [5,6]. Such an asymmetry is another consequence of 
parton energy loss and yet it is not sensitive to effects of 
initial state interactions. 
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FIG. 9. Azimuthal anisotropy in Au+Au collisions as com- 
pared to the STAR [61] 4-particle cumulant result. 

In the parton model, the azimuthal angle dependence 
of the parton energy loss will be given by Eq. (14) for 
non-central heavy- ion collisions. Using this anisotropic 
energy loss in the effective modified fragmentation func- 
tions in the pQCD parton model of high-pr hadron spec- 



From both single and di-hadron spectra and their az- 
imuthal anisotropy vi (pt ) , the extracted average energy 
loss in the parton model calculation for a 10 GeV quark 
in the expanding medium is (dE/dL)id ~ 0.85 ± 0.24 
GeV/fm, which is equivalent to dEo/dL w 13.8 ± 3.9 
GeV/fm in a static and uniform medium over a distance 
Ra = 6.5 fm at an initial time tq = 0.2 fm. Compared 
to the energy loss extracted from HERMES data on DIS, 
this value is about 30 times higher than the quark energy 
loss in cold nuclei. Since the parton energy loss in the 
thin plasma limit is proportional to the gluon number 
density, one can conclude that the initial gluon density 
reached in the central Au + Au collisions at 200 GeV 
should be about 30 times higher than the gluon den- 
sity in a cold Au nucleus. This number is consistent 
with the estimate from the measured rapidity density of 
charged hadrons [64] using the Bjorken scenario [65] and 
assuming duality between the number of initial gluons 
and final charged hadrons. Given the measured total 
transverse energy dETjdr\ w 540 GeV or about 0.8 GeV 
per charged hadron [66] in central Au + Au collisiona at 
\[s = 130 GeV, the initial energy density in the formed 
dense matter is at least 100 times higher than in cold 
nuclear matter. 

The above analyses of RHIC data on jet quenching are 
based on an assumption that the dominant mechanism of 
the jet quenching is parton energy loss before hadroniza- 
tion. It is reasonable to ask whether leading hadrons 
from the jet fragmentation could have strong interaction 
with the medium and whether hadron absorption could 
be the main cause for the observed jet quenching [67,68]. 
A detailed analysis of jet quenching data concludes [69] 
that the data are not compatible with such a scenario 
of hadronic absorption and that the observed patterns 
of jet quenching in heavy-ion collisions at RHIC are the 
consequences of parton energy loss. 
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In addition to the long hadron formation time (it could 
be about 30-70 fm/c for a 10 GeV pion) estimated from 
uncertainty principle, there are many patterns of the sup- 
pression that can rule out the scenario of hadron absorp- 
tion. Foremost among them, the large vi at high-pj 1 and 
the same-side di-hadron correlation cannot be reconciled 
with the hadron absorption scenario. 

Since the spectra azimuthal anisotropy is caused by 
the eccentricity of the dense medium, which decreases 
rapidly with time [70] due to transverse expansion, it 
has to happen at a very early time. By the time of a few 
fm/c before the formation of any hadron, the eccentricity 
is already reduced to a non-significant value. Any further 
interaction cannot produce much spectral anisotropy. 

As shown in Fig. 7, while the away-side jet is sup- 
pressed in central Au + Au collisions, the same-side di- 
hadron correlation remains almost the same as in p + p 
and d + Au collisions. This is clear evidence that jet 
hadronization takes place outside the dense medium with 
a reduced parton energy. A recent study [71] of the 
di-hadron fragmentation functions shows that the con- 
ditional (or triggered) di-hadron distribution at large z 
within a jet is quite stable against radiative evolution (or 
energy loss). On the other hand, a hadron absorption 
mechanism will suppress both the leading and secondary 
hadron in the same way and should lead to the suppres- 
sion of the same-side correlation. One may argue that 
the jet whose leading hadron is the trigger hadron could 
be emitted from the surface and suffer no energy loss and 
therefore no suppression of the secondary hadron. This 
is not true. Even though the same-side jet is produced 
close to the surface due to trigger bias, the parton jet 
on the average still loses about 2 GeV energy [69] , which 
should be carried by soft hadrons in the direction of the 
triggered hadron. If one collects all the energy from jet 
fragmentation for a fixed value of trigger hadron px, this 
energy in central Au + Au collisions should be larger than 
in p +p collisions by about 2 GeV. This has indeed been 
confirmed by preliminary data from STAR [15], as shown 
in Fig. 10, where the scalar sum of the pr on the same 
side of the triggered hadron is indeed about 2 GeV larger 
in central Au + Au than in p + p collisions. The analysis 
[15] also sees modification of the fragmentation function 
both along and opposite the trigger hadron direction. 

Finally, if hadron absorption can suppress high-p^ 
hadrons and jets, it would most likely happen in heavy- 
ion collisions at the SPS energy. Hadron spectra at this 
energy are very steep at high-p-r and are very sensitive 
to initial transverse momentum broadening and parton 
energy loss [45]. However, the measured it spectrum in 
central Pb+Pb collisions only shows the expected Cronin 
enhancement [72,73] with no sign of significant suppres- 
sion. More recent analysis of the Pb+Pb data at the SPS 
energy also shows [74] that both same-side and back-side 
jet-like correlations are not suppressed, though the back- 
side distribution is broadened. 
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FIG. 10. Charged hadron multiplicity (a) and total scalar 
Pt (b) of the near and away side jet with 4 < p^, %a < 6 GeV/c 
and 0.15 < pr < 4 GeV/c in p+p (open symbols) and Au+Au 
collisions from STAR experiment [15]. 



VI. SUMMARY AND OUTLOOK 



In summary, with the latest measurements of high-p^ 
hadron spectra and jet correlations in p + p, d + Au and 
Au+Au collisions at RHIC, the observed jet quenching in 
Au + Au collisions has been established as a consequence 
of final-state interaction between jets and the produced 
dense medium. The collective body of data, suppression 
of high-j»T spectra and back-to-back jet correlation, high- 
Pt anisotropy and centrality dependence of the observ- 
ables, points to parton energy loss as the culprit of the 
observed jet quenching. A simultaneous phenomenologi- 
cal study within a LO pQCD parton model incorporating 
the parton energy loss describes the experimental data of 
Au + Au collisions very well. The extracted average en- 
ergy loss for a 10 GeV quark in the expanding medium 
is equivalent to an energy loss in a static and uniform 
medium that is about 30 times larger than in a cold nu- 
cleus. This leads us to conclude that the initial gluon 
(energy) density is about 30 (100) times of that in a cold 
nuclear matter. It is inconceivable to imagine within our 
current understanding of QCD that any form of matter 
other than quark gluon plasma could exist at such a high 
energy density. 

A vast collection of data on bulk properties of particle 
production in heavy-ion collisions provides further evi- 
dence of QGP formation at RHIC. Most striking among 
these data is the elliptic flow measurement [75] that was 
found to saturate the hydrodynamic limit and exhibit 
all the signature behaviors of a collective hydro flow, in- 
cluding the splitting in V2(pt) of hadrons with different 
masses [18,19,76]. Together with jet quenching, they pro- 
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vide a collection of compelling evidence for the formation 
of a strongly interacting quark gluon plasma at RHIC. 
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FIG. 11. High-pT two-hadron correlation in azimuthal an- 
gle for different orientation of the trigger hadron with respect 
to the reaction [14]. 

The discovery of jet quenching and its central role 
in characterizing the strongly interacting quark gluon 
plasma will usher in another new era of studying the 
properties of the dense matter. With higher-precision 
data extending to higher pr, one can find out the miss- 
ing energy of the suppressed jet [77] and map out in detail 
the modification of the fragmentation functions and the 
energy dependence of the energy loss. The preliminary 
results from STAR [15] on modified fragmentation func- 
tions in Au + Au collisions have already demonstrated 
the potential of the study and could even shed light on 
parton thermalization and hadronization. The ultimate 
measurement of a modified FF in heavy-ion collisions will 
be the direct-7-triggered FF. Recent preliminary results 
from PHENIX [16] indicate that such a measurement will 
be within reach in the next few years. 

One can also study the di-hadron correlations with re- 
spect to the reaction plane. Shown in Fig. 11 are the 
preliminary results from STAR [14] on di-hadron cor- 
relation in azimuthal angle for different orientations of 
the trigger hadron with respect to the reaction plane in 
Au + Au collisions. This could eventually measure the 
length dependence of the parton energy loss. Measure- 
ments of charmed mesons can provide additional tests of 
the parton energy loss scenario of jet quenching, since 
recent theoretical studies [78-81] show many unique fea- 
tures of heavy quark energy loss that arc different from 
those of light quarks and gluons. 
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